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Other  Supporting  Agencies:  National  Science  Founda 
Summary  of  Research  Activities : 

(1)  Theoretical  Study  of  Ionic  Transport  in  Solid 
6-Alumina 

In  order  to  understand  the  dynamics  of  carrier  ions  (mobile  cations) 
in  6-alumina,  we  have  calculated  the  potential  energy  curves  for  correlated 
motions  of  these  ions. 

Since  the  x-rav  and  neutron  diffraction  studies  of  Na  and  Ag  6-alumina 
show  that  these  ions  move  on  a  two-dimensional  network  of  hexagons,  we 
restrict  the  motion  of  carrier  ions  to  the  lines  of  this  network. 

Our  calculation  shows  two  sets  of  alternating  vertices.  (See  Fig.  I) 

One  set  we  have  designated  as  Site  1  vertices;  these  coincide  with  the  one 
particle  potential  energy  minima.  The  other.  Site  2  vertices,  do  not 
correspond  to  the  potential  minima. 

In  order  to  calculate  the  potential  energy  surface  we  consider  (1) 

Coulomb  potential  due  to  point  charges,  (2)  Born  Mayer  repulsive  potential 
between  ions,  and  (3)  the  polarization  energy.  The  polarization  energy  is 
calculated  in  a  self-consistent  manner. 

Calculations  have  been  done  for  various  carrier  ions  (i.e.  Li,  Na, 

K,  Ag,  P,b)  with  the  assumption  that  the  crystal  parameters  remain  un¬ 
changed. 

The  composition  of  the  crystal  with  all  Sites  1  filled  with  carrier 
ions  (M+)  corresponds  to  M20‘11  A1?03  which  we  may  call  the  "ideal  structure." 
In  this  structure  the  potential  energy  of  Site  2  is  higher  than  that  of 
S I  to  1  hv  more  than  2  cV  for  all  carrier  ions  considered.  6-alumina  samples 
reported  In  the  literature  have  excess  M+  tons  and  charge  neutrality  is 
maintained  either  by  extra  oxygen  ions  or  bv  vacancies  of  aluminum  ions,  or 
bivalent  Ions  substituted  for  aluminum  ions. 
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lu  the  "Ideal  structure"  the  barrier  height  for  the  motion  ol  a 
carrier  ion  is  expected  to  be  an  order  of  magnitude  larger  than  tin* 
experimental Iv  measured  values  of  the  activation  energy  for  diffusion 
and  electrical  conduction.  This  suggests  some  kind  of  correlated 
motion  of  carrier  ions  in  the  presence  of  M+  ions  in  excess  of  the 
ideal  structure.  As  the  initial  configuration  of  ions  in  our  model, 
an  extra  M+  ion  is  placed  at  a  Site  2  in  the  ideal  structure.  As  the 
ion  at  a  Site  1  is  moved  toward  an  empty  Site  2,  the  M+  ion  at  the 
nearest  Site  2  moves  in.  (See  Fig.  2)  Potential  energy  values  are 
calculated  as  the  first  ion  is  moved  and  the  position  of  the  second 
ion  (started  at  Site  2)  is  adiusted  .o  minimize  the  total  potential 

O 

energy  of  the  system.  This  calculation  is  done  at  intervals  0.16  A 
until  the  first  ion  reaches  the  originally  empty  Site  2.  From  this 
calculation  one  can  see  that  the  minimum  potential  energy  corresponds 
to  the  two  M  ions  located  between  two  different  sites  (i.e.  between 
the  lattice  vertices).  The  potential  energy  difference  between  the 
maximum  and  the  minimum  mav  be  interpreted  as  the  activation  energy  in 
this  correlated  motion  of  the  pair. 

To  improve  the  calculation  we  considered  the  correlated  motion  of 
six  M  ions.  The  initial  configuration  is  identical  to  that  of  the 
previous  calculation.  As  one  M  at  Site  1  is  moved  toward  an  emptv 
Site  2,  the  positions  of  the  M+  ion  at  the  nearest  Site  2  and  four  M+ 
inns  at  nearbv  Site  1  are  allowed  to  adjust  to  minimize  the  potential 
energy  of  the  system.  (See  Fig.  3)  The  potential  energy  values  as  a 
function  of  the  displacement  of  the  M+  (originally  at  Site  2)  are  given 
in  Fig.  4.  Again  one  notices  that  the  potential  energy  minima  correspond 
to  the  location  of  two  ions  in  between  adjacent  vertices  of  a  hexagon. 

A  very  interesting  result  of  this  calculation  is  that  other  four  M+ 
ions  do  not  move  away  from  their  original  positions  (Sites  1),  but  each 
remains  near  its  original  location  (Site  1)  with  small  displacements 
as  shown  in  Fig.  3  while  two  M  ions  move  the  length  of  one  side  of  a 
hexagon.  The  energy  difference  between  the  maximum  and  minimum  positions 
mav  be  interpreted  as  the  activation  energy  and  these  values  are  listed 
in  Table  I  and  compared  with  reported  activation  energies  of  electrical 
conduct  Ion .  Ihe  agreement  seems  to  be  satisfactory  except  for  l.i  +  .  The 
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disagreement  on  Li+  Is  not  surprising  since  Li+  ion  is  not  expected 
to  move  along  the  sides  or  hexagons.  (A  separate  numeiical  computation 
of  a  one  particle  potential  energy  curve  along  lines  perpendicular  to 
the  hexagonal  planes  shows  that  Li  has  the  potential  energy  mimina  awav 
from  such  plane.) 

We  have  found  that  the  larger  value  of  Ag+-ion  polarizability  yields 
the  lower  activation  energy,  and  the  larger  radius  of  Na+-ion  leads  to 
the  higher  activation  energy. 

The  random  walk  model  is  used  to  calculate  oQ  for  crystals  with 
16%  excess  M  ions  and  the  results  are  compared  with  some  reported 
measured  values  in  Table  II. 

We  have  also  carried  out  similar  calculations  for  some  crystals 
with  a  low  concentration  of  Impurity  ions.  These  are  listed  in  Table  III. 
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Table  I. 

Calculated 

and  experimental 

ac.tivat  ion 

energies  of 

beta-alumina  crystals 

in  eV.  Ea(2)  and 

Ea(6)  are 

calculated  with 

2  and  6  ’l+ 

ions,  respectively,  allowed 

to  ad-just 

their  positions. 

Mate  rial 

Ea(2) 

Ea(6) 

Ea(ExP) 

Li -bet a 

0.29 

0.10 

0.378 

Na-bet a 

0.27 

0.13 

0.165 

Ag-beta 

0.15 

0.10 

0.176 

K-be  ta 

0.32 

0.23 

0.233 

Rb-beta 

0.52 

0.39 

0.311 

Cs-beta 

1.15 

0.81 

Table  II 

.  Calculated  and  experimental 

values  of  uD  for  beta-alumina 

crvstals 

with  n=0.16N.  The  values  are 

in  units  of  10  :  (ohm-cm)- 1 °K . 

Material 

0o(Cal) 

0o(ExpJ 

Li -beta 

6.43 

—  — — 

Na-bet  a 

2.33 

2.5 

Ag-beta 

1.11 

1.6 

K-be t a 

2.13 

— 

Rb-beta 

1 .66 

— 

Table  III.  Calculated 

activation  energies 

;  in 

eV  for  K+  -  and  Na+  -  be 

alumina  with  impurities.  Six  metal  ions 

are 

allowed  to  adjust  their 

positions . 

Material 

Comp  lex 

I  on 

Ea(6) 

Pure  Na-beta 

+2 

2 

Na+ 

0.13 

,  4- 

4-  4- 

+ 

IK  In  Na-beta 

Na  -K 

K 

0.22 

+ 

Na 

0.10 

Pure  K-beta 

K+-K+ 

K+ 

0.23 

4- 

.  +  + 

+ 

1  Na  in  K-bc*ta 

Na  -K 

K 

0.31 

+ 

Na 

0.10 
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Publications : 


L.  L.  Dunn,  L.  J.  LaFountain,  and  R.  E.  Jackson,  "Porosity  Dependence 
and  echanism  >f  Brittle  Fracture  in  Sandstones,’  Jour.  Geophys.  Res., 
v.  78,  n.  1L,  p.  2L03-2L17,  (1973). 

P.  J.  Poper  and  D.  E.  Dunn,  "Superposed  Deformation  and  Polymetamorohism, 
(1973)^  Z°ne’  South  Carolina"  Geol.  Soc.  America,  Bull.,  v.  8L,  p.  3373-338 

L.  J.  LaFountain  and  D.  E.  Dunn,  "Anisotropy  and  the  Coefficient  of  Sliding 
friction,  Eos,  Amer.  Geophys.  Un.,  Trans.,  v.  55,  n.  k ,  p.  L28,  (197L). 

R;  fftckson  J*  E*  "Experimental  Sliding  Friction  and  Cataclasis 

of  Foliated  Rocks,  Int.  J.  Rock  Vech.  Min.  Sci.,  v.  11,  (197L  in  press). 

L.  J.  LaFountain  and  D.  E.  Dunn,  "Effect  of  Anisotropy  on  the  Coefficient 

fri(:tlon  in  Schistose  Rocks,”  Int.  J.  Rock  Me ch.  Min.  Sci., 

(197«  :n  press). 

Summary  of  Research  Activities: 

1  Dilatancy  -  Fluid  Diffusion  Model  for  Earthquake  Prediction. 

Recently  Nur  (Geology,  May  197M  has  published  confi rmation  of  th< 

1!  latency  -  fluid  diffusion  model  based  on  data  from  the  Matsti.hiro,  Japan 
earthquake  swarm.  His  interpretation  of  the  data  depends  upon  assumptions 
regarding  the  relation  between  discharge  and  pore  pressure  in  a  dilatinr 
fault  zone.  Currently  we  have  underway  a  series  of  experiments  which  ghoul < 
confirm  or  deny  Nur’s  model.  These  consist  of  fracture  tests  in  two  modes. 
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In  the  first  ve  hold  pore  pressure  constant  so  that  water  inflow  exceeds 
the  rate  of  pore  pressure  increase  (q>>  d<*/dt);  while  in  the  second  we 
allow  pore  pressure  to  decrease  so  that  d*/dt  >>q.  Within  a  few  months  we 
should  have  firm  experimental  evidence  with  which  to  refine  possible  earth¬ 
quake  prediction  schemes. 

( -  )  Brittle  Fracture  Mechanisms. 

There  follows  the  abstract  of  the  paper  cited  above  by  Dunn, 

LaFountain,  and  Jackson. 

Porosity  Dependence  and  Mechanism  of  Brittle  Fracture  in  Sandstones. 
Brittle  fracture  tests  of  105  fine-grained  quartz  arenites  were  conducted 
at  25°C,  1.0  kilobar  confining  pressure,  a  constant  strain  rate  of  6.5  X  10 ~5/acc. 
and  pore  pressure  ranging  from  0  to  750  bars.  Orientation  of  planar  anisotrony 
(bedding  or  cross-bedding)  with  respect  to  principal  stresses  has  little 
influence  on  the  fracture  strength.  The  Donath  orientation  effect  depends  upon 
rock  type.  Stronv  dependence  of  fracture  strength  on  porosity  is  of  the  form: 
y*axb  (where  y*  stress  difference  at  failure,  x*  porosity,  and  a  >0>b;  in  our 
samples  values  for  a  ranged  between  l6  and  25  kilobars,  and  b  between  -0.8  and 
-1.0).  Through-going  shear  fractures  result  from  coalescence  of  grain  boundary 
crack3,  extension  fractures  within  grains,  and  void  space.  Rocks  with  low 
porosity  develop  through-going  shears  only  after  many  grains  are  extension 
fractured.  The  functional  relationship  between  porosity  and  fracture  strength 
derives  from  the  lower  energy  required  for  propagating  cracks  to  utilize  void 
space  rather  than  forming  extension  fractures. 

(3)  Frictional  Properties  of  Anisotropic  Rocks. 

There  follows  the  abstract  of  the  paper  cited  above  by  LaFountain  and 

Dunn . 

Effect  of  Anisotropy  on  the  Coefficient  of  Sliding  Friction  in  Schistose 

RockB.  Foliated  schict  samples  were  savcut  at  *45°  to  the  cylinder  axis  and 

tested  at  500  bar3  confining  pressure,  25°C,  and  lCT^sec*1  shortening  rate,  to 

examine  the  effect  of  anisotropy  on  the  coefficient  of  sliding  friction  (w3). 

Cylinders  were  cored  from  sample  blocks  in  a  variety  of  orientations  to  produce 

sawcut-to-foliation  (STF)  angles  ranging  from  0°  to  l80°  in  15°  increments. 

The  variation  of  us  with  foliation  orientation  is  sinusoidal.  Maximum 

values  of  us  occur  when  the  foliation  is  parallel  to  principal  stress  planes 

(f.TF  =  30°  -  1*5°  and  105°  -  120°;  t  *  0),  and  minimum  values  of  us  correspond 

to  planes  of  maximum  resolved  shear  stress  (RTF  =  75°  and  l65°,  i  =  t  ) 

max. 

? 
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Sliding-surface  damage  and  gouge  development  very  inversely  with  u  . 
Surface  damage  results  from  the  interaction  of  asperities  with  the  savcut 
surface,  and  failure  ana  cataclasis  of  asperities  produces  gouge.  Precursor 
events  visible  on  force-time  records,  and  partial  loading  tests  terminated 
before  megascopic  sliding,  indicate  that  surface  damage  and  gouge  are  gener¬ 
ated  before  megascopic  sliding  on  the  savcut.  Apparently,  high  resolved  shear 
stress  produces  clip  along,  the  foliation  prior  to  sliding  on  the  savcut. 
Foliation  slip  generates  gouge  at  the  savcut  interface;  and  the  effect  of  the 
gouge  is  to  lover  ug. 
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Other  Supporting  Agencies:  National  Science  Foundation 
Publications: 

"A  Nearest-Neighbor  Technique  Applied  to  the  Reduction  of  the  Number  of 
Standard  Solutes  Needed  to  Characterize  Gas  Chromatographic  Liquid  Phases," 
S.  R.  Lowry,  S.  Tsuge,  J.  J.  Leary  and  T.  L.  Isenhour,  accepted  by  J, 
Ohromat.  Sci. 

"Progressive  Filter  Network:  A  General  Classification  Algorithm," 

S.  R.  Lowry,  J.  C.  Marshall  and  T.  L.  Isenhour,  submitted  to  Science. 

"Bayesian  Decision  Theory  Applied  to  the  Multicategory  Classification 
of  Binary  Infrared  Spectra,"  H.  B.  Woodruff,  S.  R,  Lowry  and  T.  L.  Isenhour 
submitted  to  Anal.  Chem. 

Summary  of  Research  Activities 

(1)  Studies  of  Gas  Chromatographic  Liquid  Phases 

Gas  Chromatography  has  become,  perhaps,  the  most  frequently  applied 
separation  method  in  chemical  analysis.  This  technique  is  capable  of 
separating  and  identifying  components  of  a  great  variety  of  complex  samples. 
The  selectivity  of  liquid  phases  is  typically  characterized  on  the  basis  of 
retention  behavior,  which  is  desirable  because  retention  behavior  directly 
reflects  molecular  Interactions  between  liquid  phases  and  test  solutes 
(functional  probes)  under  actual  operating  conditions.  Keller  (1)  pointed 
out,  "it  now  seems  desirable  to  characterize  the  selectivity  of  mobile  and 
stationary  phases  In  terms  of  chromatographic  behavior,  using  a  minimum 
number  of  probe  solutes  to  elucidate  the  retention  forces  involved,"  lie 
also  suggested  that  the  selection  of  these  probes  might  be  made  by  using 
the  extended  solubility  parameters  proposed  by  Snvder  (2). 


// 
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In  a  recent  work  (3)  a  nearest  neighbor  technique  was  used  to  group 
liquid  employing  retention  data  taken  with  the  first  five  of  McReynolds ' 
solutes  (benzene,  1-butanol,  2-pentanone,  1-nlt. opropane  and  pyridine)  (4). 

Using  these  results  and  other  criteria,  the  authors  proposed  a  set  of  twelve 
preferred  phases  and  presented  a  substitution  table  listing  each  of  the 
22b  I lquld  phases  reported  by  McReynolds  with  its  nearest  preferred  phase 
•n  determined  bv  a  five-dimensional  distance  calculation.  A  major  problem 
which  remains,  however,  is  the  systematic  characterization  of  McReynolds, 
as  well  as  new  phases  to  be  reported  In  the  future.  Characterizing  these 
phases  using  the  five  McReynolds  probes  would  be  a  task  of  considerable 
magnitude.  To  reduce  the  amount  of  work  involved  the  feasibility  of  using 
fewer  than  five  of  McReynolds  probes  for  the  characterization  of  liquid 
phases  was  studied  using  a  nearest  neighbor  technique. 

Two  sets  of  three  test  probes  and  several  sets  of  four  probes  gave 
results  similar  to  those  obtained  using  all  five  of  the  test  probes 
evaluated  by  McReynolds.  The  interactive  forces  important  in  gas  chromatography 
were  related  to  those  solutes  present  in  the  best  sets  of  test  probes.  It 
is  interesting  to  note  that  these  results  obtained  using  the  nearest-neighbor 
method  agree  with  the  work  of  other  researchers  who  used  statistical  methods 
to  study  retention  data. 

(1)  Leary,  J.  J.  Justice,  J.  B.,  Tsuge,  S.,  Lowry,  S.  R.,  and  Tsenhour,  T.  L., 

J.  Chromatog.  Sci.  11,  201  (1973). 

(2)  Keller,  R.  A.,  J.  Chromatog.  Sci.  11,  49  (1973). 

(3)  Snyder,  I..  R.,  "Modern  Practice  of  Liwuid  Chromatography,"  J.  J.  Kirkland,  ed., 
Wiley- Interscience,  New  York,  1971,  p.  125. 

(4)  McReynolds,  W.  O.,  J.  Chromatog.  Sci.  8,  685  (1970). 

(2)  Progressive  Kilter  Network 

With  the  ever  Increasing  use  of  computers  in  all  areas  of  science, 
applications  of  pattern  recognition  have  become  common  in  many  diverse  fields 
(1,2,3).  The  goal  of  any  pattern  recognition  algorithm  is  to  classify  objects 
or  patterns  on  the  basis  of  certain  attributes  or  parameters.  Frequently, 
classification  is  accomplished  by  evaluating  those  parameters  with  a  discriminating 
function.  Unfortunately  even  after  considerable  effort  and  computation  time 
have  gone  into  developing  a  given  discriminant  functions,  still  more  lengthly 
calculations  (often  requiring  a  computer)  are  necessary  to  claseify  unknown 
patterns.  Often  the  scientist  who  could  benefit  most,  has  neither  the  back¬ 
ground  nor  the  equipment  to  apply  these  techniques  to  his  own  research. 
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This  work  describes  a  pattern  recognition,  algorithm  which  is 
conceptually  simple  and  can  be  applied  without  resortto  a  computer. 

The  technique  employed  is  a  sequential  filter  approach. 

inch  parameter  (or  dimension)  of  the  data  is  investigated  for  an 
upper  and  lower  threshold.  Thresholds  are  defined  as  values  which  will 
classify  only  one  class  of  data.  The  first  filter  is  selected  as  that  which 
has  tlu  two  thres  lolds  that  classify  the  most  members  of  the  training  set. 

(The  thresholds  n  ;ed  not  classify  into  the  same  categories.)  These  patterns 
are  then  removed  from  the  data  set  and  the  next  best  filter  is  sought. 

The  network  can  be  applied,  once  developed,  to  classifying  unknown 
patterns  by  following  a  simple  filtering  process  of  companing  the  magnetudes 
of  individual  parameters  to  ttie  recorded  thresholds.  Successful  applications 
were  made  to  Interpreting  mass  spectra  and  electrocardiograms. 

(1)  Isenhour,  T.  L.  and  Jurs,  P.  C.  Anal.  Chem. ,  43  (10),  20A  (1971). 

(2)  Rogers,  0.  S.  and  Tanimoto,  T.  T.,  Science  132,  1115  (1960). 

(3)  Jain,  V.  K.,  International  Journal  of  Computer  and  Info.  Sci.,  2,  231  (1973). 

(3)  Baveslan  Decision  Theory  Applied  to  the  Multicategory  Classification  of 
Binary  Infrared  Spectra 

Pattern  recognition  techniques  employed  for  the  analysis  of  infrared 
spectra  have  been  previously  reported  as  acceptable  alternatives  to  search 
and  compare  methods  (1,2).  The  advantage  of  using  pattern  recognition  is 
the  ability  to  determine  common  characteristics  of  spectra  from  similar 
compounds.  Other  compounds  displaying  similar  characteristics  in  their 
spectra  can  then  be  predicted  to  belong  to  the  same  class  of  compounds 
as  the  original  sit  in  which  these  characteristics  were  observed.  Previous 
works  in  this  field  have  reported  results  using  binary  classification 
techniques.  This  means  that  the  answers  were  simply  yes  or  no;  e.g.  yes, 
the  compound  was  a  carboxylic  acid  or  no  it  was  not.  Similar  questions 
were  asked  for  esters,  aldehydes,  and  each  of  the  other  classes  of  compounds 
being  tested.  For  a  given  spectrum,  if  after  all  the  questions  were  answered 
onlv  one  answer  was  yes,  then  a  specific  prediction  could  be  made.  Otherwise 
the  prediction  would  merely  be  that  the  compound  belong  to  one  of  the  positive 
c  Lasses . 
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Another  approach  Is  to  use  multicategory  classification  techniques. 

In  this  case,  rather  than  asking  a  number  of  yes/no  questions,  only  one 
question  is  asked.  To  which  class  does  the  compound  belong?  Selection 
of  the  proper  discriminant  -unction  results  in  a  prediction  for  the  correct 

Wt'r  l°  l"IS  queStlon-  7,118  PaPer  reports  on  an  investigation  of  multicategory 
predictions  using  binary  data.  With  the  knowledge  that  the  spectrum  belonged 
to  one  of  fourteen  possible  classes,  the  goal  was  to  select  the  pronor  cliss. 

The  recognition  results  obtained  using  a  Bayesian  approach  were 
considerably  better  than  the  IX  figure  expected  from  random  guessing.  Bayesian 
approache  aie  not  used  very  often,  usually  due  to  the  difficulty  of  obtaining 
P (A  C)*  kn0;’n  as  th*  a  Prlc,ri  probability.  However,  in  a  situation  such  as 
thC  n<  ,<  scribec1  here,  it  is  quite  eesy  to  determine  the  a  priori  probabilities 
from  the  training  set.  Given  that  p(A'c)  is  obtainable,  the  Bayesian  approach 
can  then  be  viewed  merely  as  a  formalization  of  cormon  sense  (3).  If  the  order 
1,  2.  3,  4,  ,  most  frequently  occurs  when  a  compound  belongs  to  class  ],  then 
it  is  common  sense  to  predict  a  spectrum  giving  that  same  order  is  of  an  acid. 

As  was  mentioned,  the  most  likely  next  step  for  this  type  of  approach  would  be 
to  use  very  large  training  sets  incorporating  fewer  restrictions.  Also,  the 
fourteen  classes  selected  for  this  study  would  most  likely  not  be  the  ones 
selected  for  a  thorough  investigation  of  a  large  data  set.  When  using  a  table 
generated  from  a  large  data  set,  if  the  percentage  of  recognition  still  was 
near  ^0*  then  one  would  have  a  quick  and  accurate  means  of  predicting  the  type  of 
compound  without  having  to  search  large  numbers  of  spectra.  This  study 
demonstrated  the  feasibility  of  employing  a  Bayesian  approach  for  the  classif icatio 
of  infrared  spectra,  but  varying  degrees  of  utility  should  be  found  when  using 
discriminant  functions  generated  from  other  types  of  data. 

fJ)  ^!:kii^9«,09U6",,.P-  C-  I"enh0Ur’  T-  ^  KCillcy*  <  •  ^1. 

(’!  bidden  III,  K.  W.,.!urs,  I*.  C.  ,  Appl.  Spectrosc.,  27  ,  37]  (1973). 

(  I)  Kelly,  |>.  Anal,  (hem.,  44  (1  1)  (1972). 
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(4)  Failure  Analysis  of  Lubricants 

A  preliminary  Investigation  of  engine  oil  data  was  made  with  the 
cooperation  ot  the  Naval  Oil  Analysis  Program  under  the  Naval  Air  Systems 
Command.  Hie  bodv  of  a  report  on  that  initial  investigation  follows. 

The  data  set  Investigated  consists  of  information  taken  from  fourteen 
distinct  aircraft  engines.  Samples  of  engine  oil  were  analyzed  for  trace 
any  -mts  ot  each  ol  ten  elements.  The  analyses  were  done  by  an  emission 
technique.  In  each  of  the  fourteen  cases  samples  were  periodically  with¬ 
drawn  and  analyzed  until  a  disabling  malfunction  occurred.  The  ultimate 
goal  of  this  investigation  is  to  determine  if  any  factors  involving  the 
given  concentrations  foretell  engine  malfunction.  To  insure  practicality 
the  factors  need  to  be  simple  enough  to  be  visually  or  graphically  seen. 

For  tiie  following  discussion  the  set  of  concentrations  taken  from 
one  oil  sample  are  called  a  pattern.  Also  the  set  of  patterns  for  one 
engine  are  designated  a  subset. 

Observations 

Consider  the  nature  of  the  data.  The  following  observations  may  be 

made . 

1.  Some  subsets  are  represented  by  too  few  patterns.  To  distinguish  trends, 
enough  patterns  must  be  present  to  scan  the  functioning  engine  before 

the  malfunction.  Arbitrarily  the  four  subsets  with  fewer  than  four  patterns 
are  excluded  in  data  analysis. 

2.  Engine  operation  times  are  not  incl.ded.  For  any  analysis  of  concentrations 
versus  time  a  measure  of  operation  time  is  needed.  Since  only  the  Julian 
date  Is  given  for  the  patterns,  derivative  methods  will  be  unreliable. 

3.  011  changes  appear  to  be  made  randomly  and  oil  samples  seem  to  be  with¬ 
drawn  the  same  way.  Concentration  profiles  have  gross  jumps  where  oil 
Is  sampled  close  to  the  time  of  oil  change. 

4.  The  ultimate  cause  of  engine  failure  is  not  given.  There  is  little 
re. i son  to  expect  factors  indicating  one  type  of  failure  to  be  the  same 
as  factors  for  a  different  cause  of  failure. 

.>.  There  is  questionable  accuracy  In  parts  of  the  data.  This  may  be 

Illustrated  by  observing  the  set  of  concentration  profiles  for  silicon. 

The  maximum  value  for  Si  is  52  ppm,  whereas  the  other  values  lie  between 
0  and  13  ppm. 

6.  'recision  may  also  he  a  limiting  factor  in  the  emission  output  data. 

As  a  rule  of  thumb  the  absolute  error  may  be  considered  the  maximum  of 
3  ppm  and  10%  of  the  concentration. 
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Three  methods  of  data  analysis  are  applied  to  the  data  set.  Motivation 
and  expectation  for  each  of  the  three  methods  are  outlined  below. 

1.  If  anv  one  element  predicts  malfunction,  then  something  abrupt  should 
happen  to  its  concentration  at  the  time  of  failure.  Thus  a  simple  look 
at  concentration  profiles  seems  reasonable.  The  difficulty  here  is 
sampling  technique  (as  mentioned  in  observation  three).  Little 
illustrative  should  result  from  this  look. 

2.  To  overcome  the  sampling  problem,  perhaps  one  of  the  elements  may  be 

set  as  a  standard.  This  suggests  a  ratio  method - a  look  at  concentration/ 

concentration  profiles.  It  is  hoped  that  as  an  engine  malfunctions, 
the  valt-'  of  the  ratio  reaches  a  maximum/minimum  value  and  that  the 
numeric  .’•»  1  ue  of  the  ratio  is  significant.  A  limit  to  this  approach 
is  the  larp?  relative  error.  This  will  be  especially  apparent  for  low 
concentration  species.  Ratios  between  "high"  concentration  species 
are  preferred. 

3.  The  third  procedure  assumes  that  getting  to  the  extreme  is  less 
important  than  how  the  extreme  is  approached.  Simply  this  implies  a 
derivative  or  slope  technique.  For  the  given  data  set  two  immediate 
problems  exist.  Since  tie  derivative  operation  is  especially  sensitive 
to  noise,  the  large  relative  error  eliminates  the  apparent  utility  of 
this  method.  Further  theie  is  no  appropriate  time  function  available 
for  the  ordinate. 

A  fourth  procedure  might  be  a  spectral  method.  However  the  quantity 
of  data  required  for  each  subset  makes  this  impractical. 

T  rends 

Considering  the  preceding  argument,  the  concentration  ratio  procedure 
(method  two  above)  is  preferred.  Observation  of  concentration-concentration 
ratios  are  made  for  the  elements  of  each  subset.  The  denominator  is  limited 
to  high  concentration  elements  (Cu, Fe , A1 ,Mg)  and  all  of  the  36  remaining  ratios 
are  observed.  A  positive  response  occurs  when  for  elements  X  and  Y, 


([XJ) 


maL  f  unct . 
t  i  me 


max 


[xh 

[Y  ] 


subset . 


l’«  I  n  k  this  criterion  on  the  ten  subsets,  (fu)/[Fej  nn<l  iFeJ/lAl)  give  five 
positive  responses.  |t'u]/[Al]  and  (Cu]/[Mg]  give  three  positive  responses. 

Considering  all  the  problems  Inherent  in  the  data  these  results  are 
encouraging.  Unfortunately  although  the  maxima  appear,  the  value  at  these 
extrema  do  not  seem  to  be  significant. 

At  this  point  the  analysis  pleads  for  tine  factors  and  a  curve  smoothing 
algorithm.  flits  would  allow  a  derivative  procedure  to  be  attempted. 

It  needs  to  be  emphasized  that  all  the  afore  mentioned  results  only 
indicate  possible  trends.  Until  an  investigation  of  a  larger  data  set 
taker  place  all  results  are,  at  best,  tentative. 

Recommendations 

in  order  to  determine  whether  trends  are  present  or,  if  present, 
significant,  other  (actors  need  to  be  included  in  the  data.  As  mentioned 
previously  the  engine  run  times  between  sample  withdrawals  and  the  eventual 
cause  of  failure  are  needed. 

Other  suggestions  concern,  either  directly  or  indirectly,  the  oil 
change  itself.  For  data  analysis  it  would  he  ideal  If  the  oil  were  changed 
at  regular  run  time*  intervals  and  oil  sampled  immediately  before  the  oil 
change.  It  seems  that  the  latter  of  these  Is  a  reasonable  suggestion  and 
that  the  former  Is  an  ideal  expectation  Some  of  this  problem  could  be 
aided  by  including  "standard  engines"  in  the  data.  The  "standard  engine" 
would  have  its  oil  changed  regularly  and  sampled  immediately  before  the 
oil  change.  Samples  could  also  be  taken  as  a  functioi.  of  time  n-vr  the 
life  of  the  engine  oil.  These  would  provide  normal  curves  and  standardizing 
curves  to  compare  with  curves  near  malfunction  times. 

These  suggestions  intentionally  overlook  the  irreproducibility  of  the 
"coffee  can"  sampling  technique.  This  is  an  obvious  contributor  to  data 
scatter  and  need  not  he  mentioned  more. 

Uonr ins i uns 

A  lirst  analysis  of  the  data  appears  to  be  rather  inconclusive. 

Tentative  trends  appear,  hut  noise  factors  Impair  certaintv  on  such  a  small 
mil  imcomplcti  data  set.  before  a  final  decision  may  he  made  more  data  per 
pattern  ind  more  patterns  need  to  he  investigated. 
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Summary  of  Research  Activities: 


(I)  Oxidation  ol  Coordinated  Ligands 

Ih'ta  I  led  studies  have  heel!  carried 
intramolecular  electron  transfer  from  a 


out  on  several  reactions  involving 
coordinated  ligand  to  rut  hen ium ( 1  I  1 ) 


Ru  (l»  I  py  )  ^  (N.j) 


- - >  Ku(bJpy)2(N3)2f  :K  =  0.17V 

)K»(hipy)2(CH3CN)NJf  +  i/2N2:k25°=  lxlO^sec"1 


Ru  (h  i  |>y  )  2  (N.j)  2  +  Cll  jCN 


-t 


0.  hr>V 


-e 


->  Ku  (I)  I  pv )  ;  (i'll  jl'N)  N  :l  'y  , 


Kn  0’  I  pv )  ,  (fll  ji  N )  N  j  — - 
Ku(hlpv)  ,((!!  (t  N)N  +  +  Clljl’N  >•  Kh  (b  I  pv )  ,  (C!l  jCN)  ,2 '  +  i/ :  k*'  ’  -2xl<f‘\a 


(hipv  Is  2 , 2* -b ipy r id ine ;  the  E^(  values  ire  "s.  the  saturated  sodium 
chloride  calomel  electrode  in  acetonitrile) 

Ru(bipy)^ent+ - — — >  Ku(hipy)^en^+  :i‘‘\/2  =  *  •  OV 

+ 


Ru(liipy)  en 


3+ 


-3« 


->  (hipyl^Ri 


<S4 


.11  2+ 


+  411 


In  some  cases  it  has  been  possible  to  study  the  details  of  the  intra¬ 
molecular  election  transfer  step  since  the  ruthenium( 1 1 )  complexes,  e.g., 

2+ 

Rufbipyl^en  and  Ru(bipy) n (N^) ^ ,  can  he  oxidized  rapidly  to  rutheninm(Tll) 

either  chemically  (using  Ce(IV)  or  Br0)  or  el ect rochemicnl ly .  Initial 

-  *■  2- 

studies  have  shown  that  1  ,  NCS  ,  CN  ,  and  C„0.  are  also  oxidized  when 

2  4 

bound  to  bis(2 , 2 ' -bipyridine)  rutheniumd  1 1 ) . 


(2)  Electron  Transfer  Quenching  of  Excited  States. 

2+  2+  * 

The  charge  transfer  excited  state  of  Ru(bipy)^  ,  Ru(bipy)^  ,  is 

relatively  long-lived  and  luminesces  strongly  at  room  temperature  in 

solution.  We  have  studied  the  quenching  of  the  excited  state  by  a  series 

of  quenchers  using  spectrof luor imetry  and  flash  photolysis.  The  rates  of 

3+  3+ 

quenching  of  the  excited  state  by  FeOl^O)^  and  Ru(NH^)^  (in  water) 
and  bv  CII^-ffc)yH,^0^'_^,I ^  and  other  electron  deficient  organic  oxidants 
(in  acetonitrile)  are  at  or  near  the  diffusion-controlled  limit.  The 
quenching  of  the  excited  state  occurs  by  electron  transfer  quenching  and 
the  quenching  step  is  followed  by  a  slower  dark  reaction  which  takes  the 
system  back  to  equilibrium: 

24-  ?4-  * 

Ku(bipy)  f  - )  Ku(bipy)  j 

Ku(hipy)^  ♦  Ox  -)  Ru  (l>  I  py )  4-  Red 

Itu(lilpy)  ^  +  Red  -  - - >  Ru  ( b  I  py  )  2  +  +  Ox 


(Ox  and  Red 
We  hnvt  been 
photo  lysis. 


are  the  oxidized  and  reduced  lorms  ol  the 


able  to  measure  k()  by  emission  quenching 
k  values  have  been  obtained  for  reactions 


queni her . ) 
and  k  bv  flash 
wita  AO  changes 


-(,■) 


2+ 


as  ! .ity.f  as  1.7V.  Hie  rate  studies  on  roue i  ions  |n  this  At;  iegh.ii 
show  that  the  predictions  made  l>y  Marcus  theory  for  react.  Ions  In  the 
"abnormal"  tree  energy  region  are  incorrect. 

(3)  Met a  1 -Metal  Interactions  through  bridging  ligands. 

The^j-oxo-bridsed  complexes  of  ruthenium! ill),  [ (bipy) „XRu-0-RuX- 
C’ipy),]  (X  -  Cl,  N'O^),  have  been  prepared  and  characterized.  From 
electrochemical  studies  in  acetonitrile  the  ion  [ (hipy)2XRu-0-RuX(bipy) 0J 
also  exists  as  mixed-valence  +3  (Ku( III) -Ru ( i V) )  and  +1  (Ku (1 1 1 )-Ru(ll ) 7 
ions.  The  chemical,  magnetic,  and  spectral  properties  of  the  » -oxo-br idged 
dimers  can  he  explained  using  a  molecular  orbital  scheme  based  on  a  strong 
and  chemically  significant  ruthenium-ruthenium  Interaction  through  the 
bridging  oxide  ion. 

Weak  interactions  exist  between  ruthenium  ions  in  the  complexes 
[(NH3)5Rn(I.)KuCl(bipy)2i4+/3+  where  I.  Is  a  dibasic,  N-heterocycl  ic  ligand 
contain i n^an  intact  aromatic  n-system,  e.g.,  N Qs,  NQHQn,  and 

Nfi)'C=Ji<SN-  Thc  cffects  of  metal-metal  interactions  ate  seen  In  the 
electronic  spectra  of  the  ions  and  by  thc  appearance  of  intervalence 
Iranster  absorption  bands  for  the  +4  mixed-valence  ions  [ ( Nil ^ )  Ru(lII)- 
(L)Ru(iI)Cl(bipy)2]4+.  Weak  Interactions  also  exist  between  the  iron 
sites  in  the  mixed-valence,  1 , 1 ' -polyferrocene  ions  [  (r^ll  )Fe(C  II  -C  H  )- 
le(C5||4-C51|4)Fe(C5ll5)]2+/+  [(Fc-Fc-Fc)2+/  +  )  and  [  (C  H  )  Fe  (C  11  -C  H 

F^S^-S^)F-S^-S1U>F^S,!5)J3+/2+/+((Fc-Fc-Fc-Fc)3+/2^V  The 
nixed-valence  polyferrocene  ions  can  exist  as  a  series  of  oxidation-state 
isomers  which  with  regard  to  the  site  of  oxidation.  For  example,  for  the 
Ion  (Fc-Fc-Fc)2  two  energetically  equivalent  isomers  exist— Fc+-Fc+-Fc 
and  Fc-Fc  -Fc  —and  one  energetically  inequivalent  isomer  exists— Fc+- 
Fc-Fc  .  The  existence  of  oxidation  state  isomers  influences  such  properties 
°l  11,0  m  1  xed-vn  1  mice  ions  as  their  electronic  spectra,  the  position  and 
intensity  ..I  lnterv.-itei.ee  Transfer  hands  in  the  infrared,  the  presence 
of  st.itlstie.il  el  lei  ts  in  measured  reduction  potentials  and  rates  of 
iutr.-nno  1  l.n  electron  transfer  between  different  iron  sites. 
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Summary  of  Resea  roll  Act  i  v  i  t  i  cs  : 

(1)  Chemistry  and  Techi nlogy  of  Highly  Reactive  Metal  Powders 

We  recently  disc  vereil  in  our  laboratories  a  new  method  for  preparing 
metal  powders  ol  exceptional  reactivity.  The  basic  process  consists  ol 
the  reduction  ol  a  metal  salt  in  an  ethereal  or  hydrocarbon  solvent  with 


an  alkali  me  t  a  1  under  an  Inert  atmosphere.  The  boiling  point  of  the  solvent 
must  exceed  the  melting  point  of  the  alkali  metal  used  and  the  metal  salt 
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to  be  reduced  must  be  partially  soluable  in  the  solvent  used.  Reduction 
under  these  conditions  yields  dark  black  n orders  of  exceptional  reactivity, 
'lagnesium  prepared  by  this  method  reacts  with  h romobenzene  in  five  minutes 
tt  -78  C,  producing  a  quantitative  yield  of  the  corresponding  Crignard 
reagent . 

We  have  extended  this  area  of  study  to  other  metals  and  the  studies 
which  have  been  supported  by  ARPA  include  studies  on  aluminum,  zinc, 
indium,  thallium,  and  titanium. 

Reduction  of  aluminum  chloride  in  xvlene  yields  black  aluminum  powders 
of  exceptional  reactivity.  Reaction  with  iodobenzene  at  110°C  is  complete 
in  5  minutes;  this  it  to  be  compared  with  the  literature  report  that 
ordinary  alumina  showings  react  with  neat  iodobenzene  at  100°C  but  require 
44  hours.  A  communication  on  this  has  been  published  and  a  full  paper  is 
in  preparation. 

Use  of  indium  metal  in  synthesis  has  had  very  limited  use  because  of 
the  low  reactivity  of  the  metal.  We  to  nd  that  reduction  of  indium  salts 
in  xvlene  with  potassium  yields  indium  powders  of  exceptional  reactivity. 
Reaction  of  this  powder  wi th  alkvl  and  ary]  iodides  yields  the  corresponding 
dialkvl  or  diary]  indium  monoiodide  compound  in  near  quantitative  vields. 

A  communication  on  this  work  was  recently  published.  Also,  we  have  found 
that  reaction  of  the  indium  powders  with  diorganomercurv  compounds  gives 
high  vields  of  the  corresponding  trialkyl  or  triaryl  indium  compound. 

Reduction  of  zinc  salts  with  potassium  in  tetrahvdrofuran  vields 
black  powders  which  also  have  unusual  reactivity.  Panid  reaction  with 
b romobenzene  to  vield  the  corresponding  nhenvlzinc  bromide  indicates  the 
high  reactivity  as  tills  reaction  had  not  been  possible  prior  to  our  work. 

As  with  magnesium,  wo  found  that  the  reduction  ol  the  zinc  salts  in  the 
presence  of  .alkali  salts,  in  p.articular  I, it',  yields  zinc  powders  of  ever 
higher  reactivity.  Tito  effect  of  the  presence  of  the  alkali  salt  on  the 
reactivity  of  the  black  powders  is  currently  under  study.  Additional 
studies  with  the  zinc  powders  have  centered  on  their  reaction  witli 
'‘-halo  esters  which  Is  known  as  t lie  Reformatskv  Reaction.  Use  of  the 
active  zinc  powders  allows  reactions  to  be  carried  out  at  much  lower 
temperatures  and  results  in  much  higher  yields. 
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III. il  Mum  not  il  has  horn  used  vorv  I  I  t  t  lo  in  ill  net  synthesis 
because  o t"  Its  low  reactivity.  Wo  have  found  that  reduction  of  thallium 
(I)  salts  with  potassium  in  xvlene  or  dlglvme  yields  black  powders  of 
exceptional  react ivi tv.  The  choice  of  solvent  Is  important  as  different 
chemistry  is  observed  depending  on  the  solvent  used.  The  high  reactivity 
of  the  thallium  powders  is  demonstrated  by  its  reaction  with  iodobenzene. 

Studies  with  titanium  are  onlv  of  a  Dreliminarv  nature.  However, 
reduction  of  TiCli,  in  xylene  yields  black  powders  which  do  react  with 
alkvl  Iodides.  This  has  never  been  observed  before.  Additional  studies 
are  under  ay  on  the  chemistry  of  the  Ti  powders,  and  in  particular  the 
use  of  these  powders  In  the  fixation  of  nitrogen. 

Extensive  studies  on  the  physical  properties  of  the  black  metal 
powders  .are  underway  to  try  to  determine  the  high  chemical  reactivity  as 
well  as  the  apparent  high  electron  transport  properties  of  the  black 
materials.  The  results  are  only  of  a  preliminary  nature  but  thev  suggest 
that  the  black  materials  are  a  composite  material  of  both  the  metal  and 
the  alkali  salt  generated  in  the  reduction. 


